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Abstract: When the minimal supersymmetric standard model is augmented by three 
right-handed neutrino superfields, one generically predicts that the neutrinos acquire Ma- 
jorana masses. We postulate that all supersymmetry (SUSY) breaking masses as well as 
the Majorana masses of the right-handed neutrinos are around the electroweak scale and, 
motivated by the smallness of neutrino masses, assume that the lightest supersymmetric 
particle (LSP) is an almost-pure right-handed sneutrino. We discuss the conditions under 
which this LSP is a successful dark matter candidate. In general, such an LSP has to 
be nonthermal in order not to overdose the universe, and we find the conditions under 
which this is indeed the case by comparing the Hubble expansion rate with the rates of the 
relevant thermalizing processes, including self-annihilation and co- annihilation with other 
SUSY and standard model particles. 



Keywords: 



Supersymmetry, Nonthermal Dark Matter. 



Contents 



[T]. Introduction |l| 

g. The Model § 
|2.1| Sneutrino masses || 



2.2 Neutrino masses 



2.3 Choice of parameters 



|3|. When is the Dq thermal? ^ 
|3.l| Boltzmann equation || 
3.2 Self-annihilation [7] 



3.3 Co-annihilation with SUSY particles 



3.4 Co-annihilation with SM particles ^ 

Nonthermal i>o |9| 

Conclusions and future directions 10 

|A] . Thermal relic density calculation 12 

|B|. Mixed sneutrino thermal dark matter 14 

IO. Relic from decay of thermal SUSY particles 15 



1. Introduction 

Neutrino oscillation experiments have revealed that neutrino masses are nonzero Q. A 
renormalizable extension of the standard model (SM) that incorporates neutrino masses can 
be obtained by introducing right-handed neutrinos Nr (at least two) which, in conjunction 
with the left-handed neutrinos vl of the SM allow neutrino mass terms. If the masses are 
assumed to be of the same order of magnitude as the mass differences, this mass scale is of 
the order of about 0.1 eV. Various proposals have been made for "explaining" such a tiny 
mass scale, each implying a certain magnitude for the Yukawa coupling Yn [g]. 

Supersymmetry (SUSY) |J stabilizes the electroweak scale in which the quadratic 
divergences in the Higgs sector due to SM quantum loops are canceled by contributions from 
new particles of opposite statistics. In particular, in the neutrino sector, a supersymmetric 
extension of the SM implies the addition (for each generation) of two new complex scalar 
fields - the left-handed sneutrino ul, and the right-handed sneutrino Nr - as partners of 
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the left- and right-handed neutrinos respectively. Nr and Nr are gauge singlets under the 
SM gauge group, and interact with other particles only through the Yukawa coupling. 

It has become evident in recent years through cosmological and astrophysical obser- 
vations that there is a significant amount of non-baryonic dark matter. If R-parity (R p ) 
is conserved, SUSY provides a stable dark matter candidate - the lightest supersymmetric 
particle (LSP) 1 . 

If SUSY is realized in nature, it must be broken, and since we do not know the exact 
mechanism responsible for SUSY breaking, we will parametrize this through soft SUSY 
breaking terms in an effective Lagrangian at the electroweak scale. SUSY breaking and 
electroweak symmetry breaking lead to mixing between the vl and Nr, and we denote the 
lightest such state as vq. We will, henceforth, explore scenarios where uq is the LSP. In 
this work we explore such a dark matter candidate, when all SUSY breaking masses are of 
the order of the electroweak scale. 

If the has a significant component of Dl which interacts through the electroweak 
coupling constants g and g' , the interaction strength is big enough for it to be kept in 
thermal equilibrium in the early universe. The final relic abundance of such a thermal 
species depends on when it freezes out, varying inversely as the interaction cross-section. 
This scenario has been studied in Refs. || @, ^ . Generally, a predominantly left-handed 
LSP with mass around 100 GeV annihilates too efficiently in the early universe resulting in 
too low a relic-density today, and reducing the annihilation cross-section by increasing the 
LSP mass runs into conflict with direct detection experiments. This leaves the possibility 
of a vl-Nr mixed sneutrino as a viable dark matter candidate ||]. Another possibility 
exists, as we will argue here: (i) if the vq is almost purely Nr, and (ii) if Y/v is sufficiently 
small, the interaction cross-section is too small for the v§ to ever be thermalized, and the 
above analyses do not apply. The question then arises if such a vq can be nonthermal dark 
matter, and we investigate in this work under what conditions this can be the case. 

The outline of the paper is as follows: In Sec. [2] we write down the most general R p 
conserving renormalizable superpotential (focusing on the lepton sector), and the most 
general renormalizable SUSY breaking terms, including terms that violate lepton-number. 
We write down the sneutrino mass matrix and derive the vl-Nr mixing relations. We 
impose the constraint from neutrino masses, and specify our choice of parameters. In Sec. || 
we write down the Boltzmann equation for the 9$ number-density, and obtain conditions 
that have to be met in order for the vq to be thermal. For this we include self-annihilation 
processes of the Vq, annihilation with other SUSY particles, and annihilation with SM 
particles. In Sec. [| we point out in what cases the above conditions are not met, and 
thereby identify when the vq can be nonthermal dark matter. We offer our conclusions 
in Sec. |5[ In App. |A|, we review in general the standard computation of the relic-density 
of a species that was in thermal equilibrium in the early universe, and apply it to the 
relic-density calculation of the well-studied mixed sneutrino case in App. [B]. Even if the v§ 
is never in thermal equilibrium, a Dq relic-density could result from the decays of heavier 
thermal SUSY particles; we compute this relic-density in App. [C]. 

In the case of R-parity violation the dark matter constraint can also be satisfied, e.g. by a very light 
and, thus, very long-lived gravitino LSP [H. 
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2. The Model 



To the field content of the MSSM, we add (for each generation) the right-handed neutrino 
superfield N = (Nr,N,Fn). Written as left-chiral fields, the superfields are: Q, U c , D c , 
L, E c , N c . As usual, the MSSM Higgs doublet superfields are H u and Hd- 
The most general R p conserving renormalizable superpotential is 

W = U c YuQ ■ H u — D c YdQ ■ Hd + N C Y^L ■ H u — E c YeL ■ Hd + N c ^-^-N c + fiH u ■ Hd , (2.1) 

where A ■ B denotes the antisymmetric product of the fields A and B, Y are the Yukawa 
couplings that are 3x3 matrices in generation space, and, M/v breaks lepton number. In 
general the Y's are arbitrary complex matrices, and Mn is complex symmetric. However, 
without loss of generality, we can rotate N c by a unitary matrix to make M/\r real and 
diagonal with this rotation matrix absorbed into a redefinition of Yjy. 
The SUSY breaking terms are 



C 



SUSY Br 



QL^L 



-A 



R m z u u R - d R m d d R + (-u R A u q L ■ h u + d R A d q L ■ h d + h.c. 



i\m\i L - N R m 2 N N R - e R m 2 e e R + (-N R A N £ L ■ h u + e\A e t L ■ h d + h.c. 



st 
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+ (bfih u ■ hd + h.c. 



(2.2) 



where q and b^M^ break lepton number, and h u and hd are the scalar components of H u 
and Hd respectively. Note that A can be an arbitrary complex matrix, and without loss of 
generality, m 2 is Hermitian while b^M^ and eg are complex symmetric. 

2.1 Sneutrino masses 



The sneutrino mass terms resulting from Eqs. (|2.1| ) and ( |2.2D are (the generation structure 
is suppressed) 
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(2.3) 



-H*v d Y N + 



v u An), and A 2 , = (m 2 z /2)cos2(5 is the D-term contribution. Based on the v mass matrix, 
we make the following observations: 



m 



2 RL = (— ^vjYn + v u An) results in z>£ <-> N R mixing. 



C£ breaks lepton number and results in ul <-> v% mixing pC 
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Mj\r breaks lepton number and results in vl *-* N R mixing. 
bftMN breaks lepton number and results in Nr <-> N R mixing. 



The phases present in Eq. (2.3) lead to CP violating effects which we do not explore 
in this work. We can write vl and Nr in terms of real fields z?i, V2-, N\ and N2 as 

H = (p\ + ii>2)/V2 , 

Nr = (Ni + iN 2 )/V2 . (2.4) 
In this basis, if m LL , Tn RR , ni RL , q, Y^Mn, and b^M^ are all real, Eq. (2.3) reduces to 



block diagonal form [[H], [Tl| as 
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2.5) 



m 2 RR + b N M N I 



We can diagonalize Eq. (|2.5D by performing the unitary rotation 



Ni 



where the mixing angle is given by 
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(2.6) 



2 (m 2 RL ± UaM^Ylv 



( m LL =F <*) 



(2.7) 



with the top (bottom) sign for i = 1 (i = 2). We will use the notation: Sj = sin(#^), 
cj = cos(0f) and = tan(0f ). 

We see from Eq. (|2.E| ) that the lepton-number- violating parameters q, 6 at and Mjv 
split the z>i <-> z>2 degeneracy, and also the Ni «-> iV2 degeneracy. We denote the lightest of 
the sneutrino mass eigenstates as vq, and its mass Mlsp- We will also assume that vq is the 
lightest supersymmetric particle (LSP), with R p conserved, so that it is a stable particle. 
For definiteness, we will assume that the LSP is in the upper-left 2x2 block. We denote 
the heavier sneutrino mass eigenstates collectively as vu- 

The hidden generation structure in Eq. Q2.5| ) in general implies new flavor changing 
neutral current (FCNC) contributions, which are experimentally constrained to be small 
fi~2[ l . These contributions are small for instance if the fermion and scalar mass matrices are 
suitably aligned. Moreover, an alignment at the high (GUT) scale is sufficient, since we 
expect that renormalization-group loop-induced violations of this alignment to be propor- 
tional to Yjv, which, as will be discussed shortly, we assume to be very small. 
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2.2 Neutrino masses 



The superpotential in Eq. ( |2.lD leads to the neutrino mass terms (iV c and v are 2-component 
spinors) 

Ca« = -N c v u Ynv - N C ^N C + h.c. . (2.8) 

When v u Yn <C Mjv, the lowest neutrino mass eigenvalue is given by the standard seesaw 
relation 

v 2y2 



„ v = ^. (2.9) 



2.3 Choice of parameters 



The masses of the light neutrinos, m u , are required to be less than about 0.1 eV to corre- 
spond to the mass scale inferred from neutrino oscillation experiments. The usual Type I 
seesaw results if Mjv ~ 10 14 GeV and Y/v ~ 0{1). However, this is not the only possibility; 
Mn could be at the electroweak scale [f^, ||, ||], or even at the eV scale [jl3|| . If lepton number 
is a good symmetry (i.e., if Mn = bjy = eg = 0), a Dirac neutrino results, and Yjy ~ 10 -12 
in order to have m v ~ 0.1 eV. 

We will parameterize An such that 

A N = a N Y N m e , (2.10) 

where is a dimensionless constant. We further assume a at to be order one, i.e., if the 
neutrino Yukawa couplings turn out to be very small (which is the case here), the neutrino 
j4jv also turn out to be significantly smaller than the weak scale. This assumption is 
automatically satisfied in the constrained MSSM (T^] (sometimes referred to as "minimal 
supergravity" MSSM). Other SUSY breaking mechanism that yield zero A-terms at a low- 
enough SUSY-breaking scale would also qualify. In this case, Eq. ( |2.7D then implies that 

si ~ Y N — a m ; a m = [a N -\ . (2.11) 

m e V mi J 

In this work we consider the case when all other relevant mass-scales are at the electroweak 
scale, including the neutrino Majorana mass (Afjy). In order to obtain m u ~ 0.1 eV, we 
see from Eq. ( |2.9D that we require Yjy ~ 10 -6 . We take the mass of Nr to also be around 

v u ,m e ~ 10 2 GeV, from Eqs. ( ^To|) and (|2~ll~l) , we have: 



the electroweak scale, i.e., (mL — b^Miy) ~ v (cf. Eq. (f2T5|)). For Yn ~ 10 6 and 



Y N ~ 10" b ; A N ~ a N ■ 100 keV ; s x ~ 10" b a m . (2.12) 

Note that, as long as a at is indeed order one (as we assume above), ajy is also expected to 
be of order one (see Eq. 2,11| ), and hence si is guaranteed to be very small. Hence, we will 



be interested in exploring such a predominantly right-handed LSP (Dq ~ Nr) as a dark 
matter candidate, with s± <C 1 and ci « 1. 
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3. When is the z> thermal? 

Depending on the vq interaction strength with itself and other particles, it can either 
thermalize if the interaction rate is bigger than the Hubble expansion rate, or remain 
nonthermal if the interactions are too weak. In this section we identify what conditions 
have to be consistently satisfied for 9q to be in thermal equilibrium with the SM thermal 
bath at some moment in the thermal history of the universe. 

3.1 Boltzmann equation 

The i>q number density n„ is governed by the Boltzmann equation 
d 



dt 



no = -3Hnt - {av) SA (n~ - n? Qe<? ) - {av) CA {n h) n^ - nc, oeq n^ eq ) + C r , (3.1) 



where the subscripts SA and CA on the thermally averaged cross-section (cry) denote 
self-annihilation and co-annihilation (annihilation with another species <f>) respectively, n eq 
denotes the equilibrium number density, and Cr is the contribution due to decay of heavier 
particles into Dq. 

Cr includes the contribution from all heavier SUSY particles that can decay into vq. 
In App. we identify such heavier SUSY particles, and, assuming they are in thermal 
equilibrium, compute the present relic density of the Dq by integrating the Boltzmann 
equation. It turns out that these contributions could contribute too big of a Dq relic- 
density, and have to be "dealt with" appropriately. In the remainder of this section we will 
consider the effects of the other terms in the Boltzmann equation. 

The number density will be driven to its equilibrium value provided that the inter- 
action rate is bigger than the Hubble expansion rate, i.e., if 

{av) SA np > 3H ; (av) CA n (f> > 3H , (3.2) 

for the self- and co-annihilation channels respectively. For annihilation of non-relativistic 
scalar particles {i>q here), we can make the Taylor expansion |l6| (av) = a + bv 2 + • • • ~ 
a + 6bT/M, where the first (second) term is the s-wave (p-wave) contribution. Using 
Eq. QA.2D for H, we get 

(*v) SA n^>s(lMg *-^\ ; {av) CA n^> Z^.Wg^^j , (3.3) 

as the condition for thermal equilibrium for the self- and co-annihilation channels respec- 
tively. 

Next, we check for various processes whether the inequality Eq. (3.3) is satisfied for 



T ~ Mlsp; if we find that equilibrium is not established for T ~ Mlsp, it will not be 
for T 3> Mlsp- This is because the cross-section at most goes like a oc 1/T 2 , while the 
number density goes like n tx T 3 , which implies that the left-hand-side of Eq. (|3.3|) goes 
like T while the right-hand-side goes like T 2 . Therefore, if the inequality is not satisfied at 
T ~ Mlsp it will never be satisfied for T 3> Mlsp- For this reason we only need to check 
if the condition is satisfied for T < -Mlsp- 
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Process 


Cross-section 


Limit 


Implied limit 


K) 


4 ( g 2 j_ g' 2 \ 2 
16tt ^ m& + M n ) 


si > 1U 




(b s ) 


4 (g 2 +g' 2 | g 2 +g' 2 Y M 2 

16vr \ Afi 1 Afi J M LSP 


si > 10~ 2 


a m Y N > 10~ 2 


(Cs) 


4^ \A N -^Y N f , 
16tt Mi JPS 

h 


Y^ Sl \A N - fi*Y N \f PS > lOkeV 


Y^a N a m Y N fps > 1CT 2 


(ds) 


y^t i t 

16n M'i A/ LSP 
Li 


Yjv > 10" 3 - 5 


y^v > io~ 3 - 5 



Table 1: Self-annihilation cross-sections and limits. We convert the limit in the third column to 
an implied limit shown in the last column using Eqs. flliol) and (|Hl|). 



3.2 Self-annihilation 

The dominant vq self-annihilation processes are 2 : 
(a s ) vqVq — > z^l^l via t-channel exchange of (Majorana) W s and -B. 
(6 S ) z>o^o — * ; W + W~ via t-channel exchange of Dh; I. 
(c s ) vqVq — ► V'V' (where ^ = (c, i, &)) via s-channel exchange of /i u ; /i^. 3 
(d s ) Po^o — ► z^l^l ; eLeT via t-channel exchange of ; H®. 

In Table |l| we show the cross-sections for the above processes for T < Mlsp, under 
the simplifying assumption that the mass of the exchanged particle is much bigger than T. 
Y^p is the Yukawa coupling of if). The limit shown in the table is got from the inequality 
Eq. ( |3,3| ), and if the limit is satisfied, that particular process can thermalize vq. The 
inequality is evaluated at T ~ Mlsp, with n„ ~ T 3 « M^gp, with all SUSY masses taken 
to be around 100 GeV, and for c\ « 1. In process (d s ), we include the factor (T/Mlsp) due 
to the p-wave annihilation, and ignore the s-wave contribution that is helicity-suppressed 
by the factor (m e /MLSp) 2 ~ 10 -4 (for the r lepton). In deriving the limit for process 
(c s ) we assume Mh u ~ Mh d = M^, and make use of the definition = {cln — (J>* /me)- 
For large tan/3, the bottom final-state can be important since Y^ = Yf, can be 0(1). fps 
denotes the phase-space factor which can be a severe suppression for the top final-state. 

For processes that limit parameters other than Yjy, we show in the last column of Ta- 
ble [l| implied limits on Y^, assuming Eqs. (2.10) and ( 2.1 1| ) . Among all the self-annihilation 



processes, (d s ) leads to the strongest limit Yn > 10 3 ' 5 . Therefore, for the case of interest, 
Eq. ( 2.12[ ), we infer that self-annihilation processes are not effective in thermalizing vq. 



2 The process, uquq —* iptp (where ip is a SM fermion) due to s-channel Z-boson exchange, is not present 
since the coupling of a pair of identical scalars to the Z-boson is zero. 

3 Associated with channels involving h u ,ha exchange, there are analogous channels with a longitudinal 
gauge-boson (Goldstone boson) exchange. We do not explicitly show these since their contributions are of 
the same order of magnitude as the h u ,hd contributions discussed. 
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Process 



Cross-section 



Limit 



Implied limit 



(g 2 +g' 2 ) s\c\ x T 



16tt 



' m£m lsp 



Cu H Sl > 10 



-6.5 



Cu H a m YN > 10 



-6.5 



i6tt ~WT j ps 



Qsifps > 10 



-6.5 



CsOtmYNfpS > 10" 



-6.5 



(Cc) 



je *t A l /j 
16tt IF 



-5 



Cf^Yjv > 10" 



(40 



gyjjP? i r 

16tt mT Mlsp 

Q. 



Ct fl *V > 10 



-6.5 



Q R Y N > 10 



-6.5 



(e c ) 



(fFggg l^-^l 2 f2 

16tt Mj •'-PS 
a 



CuhY^An - n*Y N \f PS > 10 keV 



Y,p(c, H a N Y N fps > W 



Table 2: Co-annihilation with SUSY - cross-sections and limits. We convert the limit in the third 



column to an implied limit shown in the last column using Eqs. (2.10) and (2.11) 



3.3 Co-annihilation with SUSY particles 

The Dq can thermalize by co-annihilation with other SUSY particles through processes such 
as (i/j denotes a SM fermion, and s denotes SUSY particles other than v): 

(a c ) vqvh — > ipip via s-channel exchange of Z-boson. 

(b c ) uqs — > eis' via t-channel exchange of W^. 

(c c ) VQt* R — > epbp via t-channel exchange of H . 

(d c ) i>Qt* R — ► epbp ; vlXl via t-channel exchange of H£ ; H®. 

(e c ) vqVh —> ipij> (where ip = (c,t,b)) via s-channel exchange of h u ; hd- 

We recall here that if a heavier species <p is in thermal equilibrium, its number density at a 
temperature T ~ Mlsp is Boltzmann suppressed (cf. Eq. (|A.6|) ) compared to the thermal 
number density of by the factor = e~( AM, */ T ), where AM^ = {M$ — Mlsp). 

In Table ^ we show cross-sections and limits that arise from the above co-annihilation 
processes with heavier SUSY particles. We denote the Boltzmann suppression of the heavier 
species <p as Q, (see above), the phase-space suppression factor due to heavy particles in the 
final state as fps, and the ip Yukawa coupling as K/,. Processes (a c ) and (d c ), like (d s ) in 
the previous subsection, has the p-wave annihilation factor (T/Mlsp), and we ignore the 
helicity-suppressed s-wave contribution. The implied limit if Eqs. ( |2.1C| ) and ( 2.11 ) hold is 
shown in the last column. We have defined, as before, on = (ajv — M* l m l)- 

As we see from Table |2[ though the processes impose strong constraints, if phase- 
space and Boltzmann suppression factors are strong enough, for the case of interest here 
(Eq. ( [2.12D ), we expect co-annihilation processes with heavier SUSY particles not to be 
effective in thermalizing Dq. We note that if tan(/3) is large, co-annihilation with bp can 
also be important, and for Yf, ~ 1 will lead to limits similar to that from processes (c c ) and 
(d c ). Also, in this case, in process (e c ) the bottom final state is important. 
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Process 


Cross-section 


Limit 


Implied limit 


(om) 


y^t-^r \a n -h*y n \ 2 f2 

16tt Ml 1 PS 


r^|^iV-M*^|/ps > lOkeV 


Y.^a N Y N fp S > 10~ 7 


(b M ) 


Y N r i, c l 1 *2 
16vr M- ^ PS 


Y^Yn/ps > 10- 7 


Y^Yjffps > 10- 7 


(cm) 


16tt M* J PS 


si/ps > 10" 6 - 5 


amY N f PS > 10" 6 - 5 


{d M ) 


9 4 ^ c 2 M LSP f 2 
16tt M| 'PS 


si/ps > 10" 6 - 5 


a m Y N f P s > 10" 6 - 5 




a^i Mi SP f2 
16tt Af 2 - J PS 


sifps > 10" 6 - 5 


a m Y N f PS > 10- 6 - 5 



Table 3: Co-annihilation with SM - cross-sections and limits. We convert the limit in the third 



column to an implied limit shown in the last column using Eqs. (2.10) and (2.11). 
3.4 Co-annihilation with SM particles 

The z>o can thermalize through interactions with thermal populations of SM particles. The 
important processes are (tp denotes a SM particle): 

(um) z-'oV'-R.l — ► ^h'Pl^r (where ip = (c,t,b)) via t-channel h u ; hd exchange. 
(Pm) vqiPl — ► VLtpR (where tp = (c,t,b)) via t-channel H exchange. 
(cm) voifj — > e~Lip' via t-channel exchange. 
((1m) v§tp — > vptp via t-channel Z exchange. 
(cm) ~~ > Lip via t-channel W exchange. 

Table || shows the cross-sections and the limits due to co-annihilation of uq with SM 
particles. As before, in process (om) we define = (ajv — ^* / m e), and we compute the 
limit assuming M^ u ~ M/j d = M^. In (cm), the B channel can also contribute but is less 
important given the smaller g' . 

In (ajvf) and (6a/), the charm and bottom channels are suppressed by the smaller 
Yukawa couplings compared to the top, and for T < 100 GeV the top contribution is 
Boltzmann suppressed. In general, processes (c&m) and (cm)-(^m) are suppressed further 
due to fps <C 1 if the left-handed SUSY particles in the final-state are all much heavier. 
If tan/? is large, co-annihilation with bottom, especially (&a/)> will lead to the strongest 
constraint since it is not Boltzmann suppressed (unless T < nib). However, (6a/) may be 
suppressed if Mg S> 100 GeV. Thus a combination of suppressions from fps, Cip and ^ 
can render these processes incapable of thermalizing Dq. 

4. Nonthermal z> 

If the Dq is never in thermal equilibrium, the standard relic density calculation presented 
in App. [A] cannot be applied. We consider in this section such a situation when an almost 
purely right-handed LSP (vq ~ Np) is not in thermal equilibrium, and the possibility of it 
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being nonthermal dark matter. In order for such a situation to be realized, it should be 
true that all the i>o thermalization conditions that we discussed in Sec. || are not satisfied. 
For the case specified in Eq. fl2,12| ), this happens for instance in the inflationary paradigm 
for a low reheat temperature (Trh <C 100 GeV), whereby the reheating produces only 
Nr (and SM particles), and the number-densities of heavier SUSY particles (and top) are 
severely Boltzmann suppressed. We will elaborate on this in the following. 

First, as we found in Sec. [T^, none of the self-annihilation processes can keep the Dq 
in thermal equilibrium for the parameter values in Eq. (2.12). 



Second, from Sec. 3.3, the co-annihilation processes can be important in thermalizing 



the vq for the parameter values in Eq. ( |2 .12 ). However, (b c ) is suppressed further due 



to fps <C 1 which we expect given the heavy particles in the final state, (d c ) is p-wave 
suppressed for Trh *C Mlsp, and, (a c ) and (e c ) are Boltzmann suppressed due to the 
heavy i>n\ hence we expect these processes also to be ineffective in thermalizing uq. 

Third, the processes of Sec. [l| can all be important. However, phase-space suppression 
f PS <c 1 due to the heavy particles in the final state, the Boltzmann suppression due to 
the heavy top (for Trh < M t ), and suppression due to Myh H > 100 GeV can render these 
processes ineffective in thermalizing the vq. 

Fourth, for the parameter values in Eq. ( [2.12 ), the uq relic-density from decays of 



heavier SUSY thermal particles are problematic and tend to overdose the universe if 
thermal populations of heavier SUSY particles are present (cf. App. |C|). However, for 
Trh *C 100 GeV, with no significant number-densities of heavier SUSY particles, this 
problem is avoided. 

If such a situation of a nonthermal Nr is realized, the relic abundance is determined by 
how the Nr is produced in the early universe. For example, if the universe went through 
an inflationary epoch driven by a scalar field 3> (the inflaton), the relic abundance today 
becomes directly related to the coupling of to the Nr (and the other particles), and the 
details of reheating. 

5. Conclusions and future directions 

We analyzed the sneutrino sector after writing down the most general R p conserving renor- 
malizable effective theory, including lepton-number violating terms. We considered the case 
when the neutrino Majorana mass term -Mat, and all SUSY breaking masses are at the elec- 
troweak scale, which implied that the neutrino Yukawa coupling Yn ~ 10 -6 . Assuming 
that the A-term is proportional to Y^ we were led to a tiny sneutrino mixing angle s ~ Y/v- 
We analyzed the cosmological implications of having such a predominantly right-handed 
sneutrino LSP (pQ w Nr), for values given in Eq. ( 2.12j ). 



Since such a vq interacts mostly through the tiny Yjsj, its self- and co-annihilation 
cross-sections are tiny. Therefore, the Dq should not be in thermal equilibrium at any time 
in the early universe, for if it is, the standard thermal relic-density calculation (reviewed in 
Apps. [A| and [B|) indicates that it would severely overdose the universe since it freezes-out 
too early. Moreover, we argue in App. |C] that decays into vq from heavier thermal SUSY 
particles, for these parameters, also overdoses the universe (see also fl4j| ). 
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The main focus of the paper is in investigating in which case the vq remains nonthermal, 
so that the problems mentioned in the previous paragraph for a thermal relic do not apply. 
In order to answer this, we start with the Boltzmann equation for the Dq number-density, 
and compare its self- and co-annihilations with the Hubble expansion rate; as shown in 
Eq. fl3.3p . For a given process, if the interaction cross-section is smaller than the Hubble 
rate, this process cannot thermalize the vq. We check to see whether dominant interaction 
processes of the Dq with itself, other SUSY particles, and SM particles, satisfy this condition. 

In Sec. 3.2 we find from Table [l] that the self-annihilation processes are ineffective in 
thermalizing such a Dq for the values in Eq. ( 2.12| ). Next, in Sec. 3.3, Table |2| we find 
that Dq co-annihilation processes with other SUSY particles (if they are present in thermal 
number-densities) are quite effective in thermalizing the Dq. However including phase- 
space and Boltzmann suppression factors shown in the table can render these processes 
ineffective. Lastly, in Sec. U, Table | we find that co-annihilation processes with SM 
particles are significant in thermalizing the Dq; however, the phase-space, Boltzmann, and 
heavy intermediate particle mass suppression factors can render them ineffective. 

In Sec. U we argued that in the inflationary paradigm, if the reheat temperature Trh *C 
100 GeV, all of these constraints are rendered harmless and the Nr is nonthermal. The 
relic number density will then be connected to details of reheating and the z>o coupling 
to the inflaton. Exploring these details and of how baryogenesis works in this scenario is, 
however, beyond the scope of this work 

In colliders, any heavy SUSY particle that is produced undergoes cascade decays into 
the vq LSP, which exits the detector as missing energy. An important signature associated 
with a predominantly Nr LSP is the occurrence of a displaced vertex in the detector in the 
case of Majorana neutrinos, or even a very long-lived charged particle leaving the detector 
in the case of Dirac neutrinos as the life-time scales like |Y/v|~ 2 - Further details will be 
discussed in an up-coming paper . 

In summary, if all new physics scales, including the SUSY breaking and the lepton num- 
ber breaking scales are small, there is a good chance that the LSP is the lightest sneutrino 
Dq. Furthermore, depending on the mechanism of SUSY breaking, the could naturally 
be an almost-pure right-handed sneutrino. Here, we argue that, if the vq is stable, such a 
scenario is not only cosmologically allowed, but can also properly fit our understanding of 
dark matter, as long as very stringent requirements on "initial conditions" are met. Such 
a scenario may also address some of the problems with the standard cold dark matter 
(CDM) paradigm |2(|. Moreover, collider signatures and signatures associated with such a 
dark matter candidate seem unusual enough to render such a scenario phenomenologically 
intriguing. 
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A. Thermal relic density calculation 



Observations indicate at a high confidence level that the universe is flat, and we therefore 
assume that such is the case in presenting all the formulas here. In order to derive the 
relic abundance of the dark matter candidate x> we start with the Friedmann equation for 
a flat universe 



8vrG 



-P , 



(A.l) 



which can be written during the radiation-dominated era in terms of the temperature as 

T 2 



H(T) = 1.66 5 



0.5 



Mpi ' 



(A.2) 



where g* is the effective number of relativistic degrees of freedom at T. The species freezes 
out at the temperature Tf when the self-annihilation rate of x becomes roughly equal to 
the Hubble expansion rate H, i.e., 



(av) n x ~ H{T f ) 



(A.3) 



where (av) is the thermally averaged cross-section and n is the number density of x- The 
number density is given as JL5|] 



n(T) 



(((3)/ir 2 )gT 3 (Boson), T>M, 

(3/4)(C(3)/7T 2 ) 5 T 3 (Fermion), T > M , 



(A.4) 



where g is the number of degrees of freedom of x (f° r example, g = 1 for a real scalar, 
g = 2 for a Dirac fermion, etc.). 

If x can co-annihilate with another species 4> then we can generalize Eq. ( A.3 ) to 

(av) SA n x + (crv) CA ~ H(T f ) , (A.5) 

where the subscripts SA and CA on the cross-sections denote self- and co-annihilation 
respectively. When x an d 4> are in thermal equilibrium, and furthermore if > M x , 
Eq. ( A.4 ) allows us to write 



n<t, 



Mx) 

with AM^ = (M^ -M x ). Then Eq. ( [Op becomes 



M x , 



(A.6) 



(<™)SA + ( av )cA ( JfJ 



AM^/T 



1U 



H(T f ) 



(A.7) 



when x freezes-out. 



The entropy density is given as |15|] 

2tt 2 



s(T) 



45 



g*T 3 = 0.44^T 3 , 



(Ai 
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where g* is the number of relativistic degrees of freedom at T. Defining the ratio of n to s 



n 



Y = -, (A.9) 

s 

we find, by writing the Boltzmann equation Eq. ( |3.l[) in terms of Y, that it is conserved 
(i.e. dY/dt = 0) provided: (a) the (av) term is negligible, and (b) Cr is negligible 4 . If the 
above two conditions are satisfied, the conserved value after freeze-out is given by 

nt 3.77 xt 1 . , . 

Y f = —L = I (A 10) 

; s f g* f MMpi (av) ' y ' J 

where M is the mass of \ and xt = M/Tf. The entropy density today sq is given by 

2vr 2 

So = — c/* r 3 Rs (3 x 10" 4 eV) 3 = 3678 cm" 3 , (A.ll) 

where we have used = 6.3 and To = 2.7 K = 2 x 10 -4 eV is the photon temperature 
today, from which we get the number density of n today 

n = Y f s . (A.12) 

We can then compute the present relic energy density of the species x from 



n M _ 10- 10 x/eV 3 
glfMpi Pc (av) 



where p c = (2.95 x 10 3 Vh eV) 4 is the critical energy density. With g*f m 100, we get (cf. 
Eq. (3.4) of Ref. @) 

n h 2 = W- 29 x f ( ^) = 1.1 x 10- 28 x/ ( . (A.14) 

For a species that is non-relativistic at freeze-out, i.e., cold dark matter (CDM), the 
number density is given by (cf. Eq. ( A. 4 )) 

' M TA 3/2 



n f = g 



2tt 



L e~ M / T f . (A.15) 



From Eq. (A.3), Xf is implicitly given by (cf. Eq. (2) of Ref. |T 



/ OmSgMMpi (av) \ 
x f= ln [ -oJZ=UI • (A. 16) 

V 9 *f x f J 



Xf depends logarithmically on M, as can be seen from Eq. ( A.16| ). In order to obtain the 
observed f^o ~ 0.3 (with h 2 0.5), for M 100 GeV, g = 1, we therefore need Xf w 21, 
or equivalently, (av) 1.5 x 10~ 9 GeV -2 0.45 pbarn. 



4 If Cr is not negligible see Eq. (|C.2|). 
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B. Mixed sneutrino thermal dark matter 



We found in Sec. ^ that if si > 1CP 3 then the will be kept in thermal equilibrium by 
self- and co-annihilations. A thermal relic results if this condition is satisfied, and in this 
section, we compute the present relic density resulting from such a thermal population of 
vq. This has been explored earlier in [fjj. 

With the lepton-number- violating parameters nonzero the vqVqZ coupling is zero (6|. 
This prevents the too efficient uq self-annihilation via s-channel Z-boson exchange and 
therefore leads to an adequate relic-density as we will show here. Also, due to the s\ sup- 
pression in the itr nu cleon interaction cross-section, direct-detection constraints are easily 
evaded jjj. Given that there are three generations, the uq are really three states, and 
the presence of these "co-LSP" states may lead to interesting effects depending on mass- 
splittings between them. For instance, the vqv'qZ coupling may now be non-zero (where 
the v'q is a heavier co-LSP state), and have implications in direct-detection experiments. 
We will not explore this further in this work. 

In App. |A| we gave the details of estimating the present relic abundance, and Eq. ( |A. 14| ) 
is of particular relevance here. In principle, {av) is the sum of the cross-sections of all the 
processes that we identified in Sec. ^2 and 3J3. However, for simplicity we consider the 
case when H,vhA are all much heavier in which case the only processes that are relevant 
are (a s ) and (c s ). We thus have (av) ag^ + 4a'' anc ^ f° r x f ~ the present relic 
abundance of Dq is 



n n h 2 



100 GeV\ /2 
1 +9 



M 



w 



100 GeV \l 2 Y_l ( \A N -fi*Y N \ \ 2 / 100 GeV \ 2 
Mg )\ ~M~h ) V M h ) 

(B.l) 

If Dq is predominantly right-handed (si -C 1, c\ « 1, uq k, Nr), if \An — n*Y^\ > Mh, 
and if AM$ H = {M^ H — Mlsp) is n °t too large compared to T, co-annihilation of the Nr 
with vjj by process (e c ) can be effective, with the cross-section given in Table |2| In this 
case, the present relic-density is 



n o/l 2 = L e (AM % /T) 

c l 



100 GeV 



A N -fi*Y N \ 



(B.2) 



We can see from Eq. ( ^7z| ) that in order to have si <C 1 when \An — M*ljv| ^ M^ requires 
mi 3> v . This implies that AM„ H is necessarily large, leading to too big a relic-density, 
and excluding this possibility. 

In general we find that if a thermal population of Nr is generated due to the processes 
of Sees. p.2| and 3.3, then the resulting relic density is too big. The reason for this is that 



the co-annihilation cross-section is big enough to thermalize the Nr, but is too small to 
keep it in thermal equilibrium long enough to deplete the number-density to acceptable 
levels before freeze-out. Therefore, if the vq is in thermal equilibrium we must have s\ not 
too small, as dictated by Eq. (|B.1|). 
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C. Relic from decay of thermal SUSY particles 

At some point in the history of the universe, if there exists a thermal population of heavier 
SUSY particles, they could decay into the vq giving rise to a relic abundance [14]. Here we 
focus on the case s\ <C 1 i.e., ~ Nr. The important decay channels are: 

(a£>) H u — > DqL (where L is the lepton doublet). 

Q>d) vh — > voipip via h u , (where tp is a SM fermion). 

(cd) £ — > i>oipip' via W ± . 

We find the cumulative relic abundance of Nr today by integrating the Boltzmann 
equation, Eq. (|3.1[) . We take as an estimate for Cr, 

C v ~n x T{ X ^D X) , (C.l) 

where V is the decay rate of the SUSY particle \ i n processes (od)-(cd), and X denotes 
SM particles in these processes. In this estimate, we ignore the angular dependence in the 
decays. The Boltzmann equation written in terms of Y„ is 

at s 



in the limit where the (av) term can be neglected in Eq. Q3.1| ). In the radiation dominated 
era we can change the independent variable from t to T using 

dt = ~§fi (C3) 



after which Eq. ( C.2 ) can be written as 



— Y = — (C.4) 
dT HTs y ' 



Integrating this from down to T, we have 



l'(T) - - / ' dT-^- . (C.5) 



I T HTs 

Using Eqs. ( |C.1| ), ( |A.2| ), (|A.4Q and ( |A.8| ) in the above equation we get 

TMpi f°° , e~ l l x 



J t /m x x 9 / 2 

Owing to the exponential suppression in the integrand for T < M x , Y freezes in at Yf ~ 
Y(M X ). 

The decay rate from process (an) (in the H rest-frame) is given by 

r 2 Y 2 (l-M 2 /M 2 ) 2 
r (-) ~ i^ M 6fps( M ^,M & ) , fpg(M x , My) = l (l ~ , (C7) 
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where we have taken the lepton in the final state as massless, c\ = cosO", and fps is the 
phase-space function. Due to this decay, from Eq. ( |C.6[) , we get 



CO ~ ^i Y Nf'ps^r I . (c.8) 

1V1 H JT/Mjj x ' 



Using Eqs. ( A.ll| ), ( [A, 12 ) and ( A.13| ) we find the present relic density to be 



^ {aD) h 2 -W 26 cM^f PS . (C.9) 



The decay rate of process (bo) is given by 

T ™ 256^ Mj fsPS ' (C - 10) 

where f3ps is the 3-body phase-space, and following the same procedure as in the previous 
case, we find 



2,2^2 \A N ~ H*Y N \ 2 M 9h M LSP 2 



n HbD) h' = i^( c ( - sf) ^ " M f H Hps • (en) 

The decay rate of process (cd) is given by 

Sl9 1 Hps , (C12) 



{cd) ~ 256^ 3 



and we find 



M 3 M T sp 

fi 0(c ^ 2 = 10 24 ^ 4 ^^/fp5- (C13) 

From these estimates, we infer that for Qoh 2 < 0(1), with all masses about 100 GeV, 
we need: Yjy < 10~ 13 , \Ajy — /U*Yv| ^ 10 eV, and s\ < 10~ 12 (where we ignore possible 
phase-space suppression factors). This is realized for instance for a purely Dirac neutrino, 



and has recently been considered in Ref. [14]. In this case the heavier SUSY particles are 



very long lived, and one has to worry if it decays during big-bang nucleosynthesis (BBN) 
causing photo-dissociation of the light elements, and spoiling its successful predictions. We 
find from Eq. ( |C.3j ) that in order for heavy particles to have decayed away before BBN, 
we need the life-time of any such particle to be r < Is, which implies that the width 
r > 10~ 25 GeV. Considering for example the decay of H, we find from Eq. ( |C7] ) that 
we need Yjy > 10~ 13 in order for the H to have decayed away before BBN. It is worth 



mentioning that a scenario similar to the one explored in 14] might be realizable even if 
the neutrinos are Majorana fermions. It seems, for example, that an "eV-seesaw" [13], 
which imposes ~ 1 eV and Yv ~ 10~ n might satisfy the conditions discussed above. 
A more detailed analysis of this possibility is beyond the ambitions of this paper. 
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